Monte Carlo Evaluation of Non-Abelian Statistics 
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We develop a general framework to (numerically) study adiabatic braiding of quasiholes in frac- 
tional quantum Hall systems. Specifically, we investigate the Moore- Read (MR) state at = 1/2 
filling factor, a known candidate for non-Abelian statistics, which appears to actually occur in na- 
ture. The non-Abelian statistics of MR quasiholes is demonstrated explicitly for the first time, 
confirming the results predicted by conformal field theories. 
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The quantum statistics of a system of identical parti- 
cles describe the effect of adiabatic particle interchange 
on the many-body wave function. All fundamental par- 
ticles belong to one of two classes: those that have their 
wave function unaffected by particle interchange (bosons) 
and those whose wave function gets a minus sign under 
permutation (fermions). In two dimensions, it is known 
that a number of exotic types of statistics can exist for 
particle-like collective excitations. For example, elemen- 
tary excitations of the Laughlin fractional quantum Hall 
(FQH) states exhibit "fractional" statistics: The phase 
of the wave function is rotated by an odd fraction of 
TT when two Laughlin quasiparticles (or quasiholes) are 
interchanged 2]. Even more exotic statistics can ex- 
ist when a system with several excitations fixed at given 
positions is degenerate In such a case, adiabatic in- 
terchange (braiding) of excitations can nontrivially ro- 
tate the wave function within the degenerate space. In 
general, these braiding operations need not commute, 
hence the statistics are termed "non-Abelian" . Remark- 
ably, the Moore-Read (MR) state, a state which is com- 
monly believed 0] to describe observed FQH plateaus at 

= 5/2 and 7/2 (which correspond respectively to half 
filling of electrons or holes in the first excited Landau 
level), is thought to have such non-Abelian elementary 
excitations Other possible physical realizations of 
non-Abelian statistics have also been proposed . States 
of this type have been suggested to be attractive for quan- 
tum computation [^. 

In Ref. , in order to establish the nature of the statis- 
tics of the Laughlin quasiholes, a Berry's phase calcula- 
tion was performed that explicitly kept track of the wave- 
function phase as one quasihole was transported around 
the other. Although approximations were involved in this 
calculation, it nonetheless established quite convincingly 
the fractional nature of the statistics. Unfortunately, it 
has not been possible to generalize this calculation to 
explicitly investigate statistics of the MR quasiholes 0. 
Although there has been much study of the statistics of 
the MR quasiholes in the framework of conformal field 
theories (CFT), it would be desirable toperform a direct 
calculation analogous to that of Ref. |2|. 



of this paper is to provide such a calculation, albeit nu- 
merically. Furthermore, the approach developed here is 
readily applicable to other FQH systems which are not 
easily accessible to analytic investigations. 

The evolution operator of a many-body system de- 
scribed by a Hamiltonian H{X) is in principle determined 
by the Schrodinger equation. In general, H(X) itself can 
change in time through dependence on some varying pa- 
rameter X{t). In such a case, let us define (pi{t) at a 
given time t to be an orthonormal basis for a particular 
degenerate subspace, requiring that this basis is locally 
smooth as a function of t. If A is varied adiabatically (and 
so long as the subspace does not cross any other states), 
then the time-evolution operator maps an orthonormal 
basis of the subspace at one t onto an orthonormal ba- 
sis at another t. A solution of the Schrodinger equation, 
-ipiit) = Uij{t)ipj{t), is simply given by 



(1) 



Since the matrix A is anti-Hermitian, U{t) is guaranteed 
to be unitary if its initial value U{0) is unitary. Note 
that if we vary A so that the Hamiltonian returns to its 
initial value at time t, i.e., H{X{t)) = H{X{0)), the corre- 
sponding transformation of the degenerate subspace can 
be nontrivial, i.e., ipiit) V'i(O) H- 

We explicitly demonstrate that this is the case for the 
MR state with at least four quasiholes. The analysis 
is done in spherical geometry [8j: N electrons are posi- 
tioned on a sphere of unit radius, with their coordinates 
given by (wi, wi), . . . , {upf, vn), using the spinor notation 
(i.e., u = e'-^/^ cos6'/2 and v = e'"^^"^ sin6'/2 in terms of 
the usual spherical coordinates). A monopole of charge 
25 = 2A^ -I- n — 3 in units of the flux quanta $0 = hc/e 
is placed in the center of the sphere, giving rise to 2n 
quasiholes which are put at (ui,vi), . . . , (u2n,'52n)- Us- 
ing gauge A = ($o*S'/27r)0cot 0, the MR wave function 
[3| is then given by 



^Pf = PfA, 



(a,6,...)(a,/3,...) 



The purpose where PfA, 



{a,b....){a,(3.,...) 



Hi 
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is the Pfaffian ^ of the iV x iV 
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antisymmetric matrix |9( 

[{UiVa - ViUa)(UjVa - VjUa) X 
{UiVb - ViUb){UjVi3 - VjUp) X ••• + (i <-> j)] ■ 

Pfaffian wave functions l(2Jl were first constructed in 
Ref. as CFT conformal blocks. This MR state is 
the exa ct g round state for a special three-body Hamil- 
tonian [li| and is also thought to pertain for realistic 
two-body interactions in the first excited Landau level 
Q . The presence of quasiholes in the ground state is dic- 
tated by the incommensuration of the flux with the elec- 
tron number. Physically, the MR state can be thought 
of as p-wave BCS pairing of composite fermions (CF's) 
at zero net field with quasiholes being the vortex ex- 
citations 13, Ull H^ . Each quasihole has charge e/4 
and corresponds to half a quantum of flux (because of 
the paired order parameter jSj). Eq. ((SJ describes a 
state with quasiholes created in two equal-size groups: 
{ua,Va),{ub,Vb),. . . and {ua,Va),{ui3,vp), . . .. Different 
quasihole groupings realize a space with degeneracy 2"^^ 
|l3l Il4l |. (Even though there are 2n!/2(n!)^ ways to ar- 
range 2n quasiholes into 2 groups of n, the resulting wave 
functions are not all linearly independent.) In the pres- 
ence of finite-range interactions, the exact degeneracy 
may be split by an amount exponentially small in the 
large vortex separation this case, infinitely slow 

braiding will not exhibit non-Abelian statistics, although 
for a very wide range of intermediate time scales, such 
statistics should apply 0|. The effects of disorder on 
the statistics are only partially understood 

Consider an orthonormal basis ^pi, with i — 
1, . . . , 2"~^, for the subspace with 2n quasiholes, which is 
locally smooth when parameterized by the quasihole co- 
ordinates. In order to determine the braiding statistics, 
we find the transformation Lpi Uijipj under the evolu- 
tion operator after two of the quasiholes are interchanged 
while the others are held flxed. The unitary matrix Uij is 
obtained by first solving Eq. ^ and then projecting the 
final basis onto the initial one. (Since we require ipi to 
be only locally smooth, the basis itself can nontrivially 
rotate after the quasiholes return to their original posi- 
tions). Eq. is integrated numerically: The differential 
equation is discretized and the wave-function overlaps 
(the right-hand side of the equation) are evaluated using 
the Metropolis Monte Carlo method. The computational 
errors are easily evaluated by varying the number of oper- 
ations. We aim the calculation at addressing the follow- 
ing questions: (1) What is the Berry's phase accumulated 
upon quasihole interchange due to the enclosed magnetic 
flux and due to the relative statistics? (2) What is the 
transformation matrix for the ground-state subspace cor- 
responding to the braiding operations? In the following, 
we will first describe the numerical method, then present 
the results, and compare them to CFT predictions [3.ll3|. 



In order to integrate Eq. numerically, the quasihole 
interchange is performed in a finite number of steps. If 
t/'^'-' is the value of the transformation matrix at the Ith 
step, then at the next step 

jjii+i) = [/(0[i + ^(0/2][i_ a('V2]~\ (3) 

where a[^^ = + V^fVl^^^ " ¥'?>/2. Our choice 

of the finite-clement scheme |3l will become clear later. 
In practice, in general we do not know an orthonor- 
mal basis for the MR states (01 in an analytic form, 
but we can numerically orthonormalize a set of 2"~^ 
linearly- independent Pfaffian wave functions ippu. Let 
^ij^ — [V'pfi: V'pfj] denote the normalized overlaps of dif- 
ferent states. (It is implied here and throughout the pa- 

per that [<Upy] ^ is evalu- 

ated numerically.) We then easily show that 

^(0 ^ [y(0]tv[/(0-^(;+i) /2 _ H.c. , (4) 

where w'^''^ = [V'pk, V'pf^"^^] and V^'-^ is defined by 
[^(Ojt^COp^CO = 1^ constructing an orthonormal basis 
ip-P = VjPippjj. We require F^'^ to be locally smooth as 
a function of the quasihole coordinates: The basis can 
continuously transform while the quasiholes are moved, 
but, e.g., sudden sign flips are not allowed. 

According to Eq. A^^^^ is anti-Hermitian, so that 
the transformation is guaranteed to be unitary if 

L/''-' is unitary. This explains our choice lj2Jl for discretiz- 
ing Eq. . Another feature preserved by our numerical 
scheme is that making a step forward, ippj^ — > ''Pp^^\ fol- 
lowed by a step backward, V'pft^'' ~^ "^pw results in a triv- 
ial transformation. We start at U'^^^ — 1 and flnd 
after performing ris + 1 steps for braiding of two quasi- 
holes (jis is increased to convergence). Because f/'pfi'' is 
some nontrivial linear combination of "^pfii finally, 
have to project the transformation onto the initial ba- 
sis: ^ C/("=)Cl^, where O = and 
^ij = [V'pfi J V'pfJ "* ] ■ The resulting unitary transformation 
matrix U then gives a representation of the braid group 
for quasihole interchanges. In the following, we describe 
our numerical experiments. 

The space describing 2n = 2 MR quasiholes is nonde- 
generate, so non-Abelian statistics cannot occur. There 
is, nevertheless, a Berry's phase accumulated from wrap- 
ping these quasiholes around each other. Our calcula- 
tion of this phase for the MR state is analogous to the 
one performed in Ref. 0] for the Laughlin state, except 
that our calculation is numerical and therefore requires 
no mean-field approximation. Let us first briefly recall 
results for the Laughlin wave function at filling factor 
V — 1/p. In the disk geometry, the Berry's phase x cor- 
responding to taking a single quasihole around a loop is 
given by 2it for each enclosed electron, i.e., x = 27r(A^), 
where (iV) is the expectation number of enclosed elec- 
trons j3|. Therefore, when another quasihole is moved 
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FIG. 1: Berry's phase x for looping one MR quasihole around 
the equator with another quasihole fixed at a zenith angle ^. 
A'^ = 4, 8, 16, 32, 64 is the number of electrons. The dashed 
line, x/'"" — ~l/8, shows a naive prediction. For cos^ ~ 0, 
the two quasiholes approach each other very closely and we 
see strong finite-size oscillations in the Berry's phase. For 
larger A'^ and cos^ (i.e., larger quasihole separation in units 
of the magnetic length), x appears to be converging toward 
zero. x(~cos^) = — x(cos^). 



inside the loop, the phase x drops by lix j-p which implies 
fractional statistics of the quasiholes. In spherical geom- 
etry , the same result holds unless the south and north 
poles (which have singularities in our choice of gauge) 
are located on different sides of the loop. In the latter 
case, the Berry's phase is given by x = T^{Nin — No^^t), 
where iVin(out) is the number of electron inside (outside) 
the loop. If a single Laughlin quasihole is then looped 
around the equator, its Berry's phase vanishes, but if 
another quasihole is placed above or below, the phase 
becomes x = ±7r/p. We check our Monte Carlo method 
by reproducing these results numerically. The charge of 
the MR {i/ = 1/2) quasihole is e/4, so that by analogy 
with the Laughlin state one might naively expect that 
the Berry's phase for looping one quasihole around the 
equator with another fixed above or below it is given by 
X = ±7r/8 ^3 (with an extra factor of 1/2 due to MR 
quasiholes corresponding to only half of the flux quan- 
tum). In Fig. n we show numerical calculation of x for ^ 
MR system having 2 quasiholes, one looped around the 
equator and the other held fixed. If the two quasiholes 
approach each other too closely, we see strong finite-size 
oscillations in the Berry's phase. However, for larger sep- 
aration, X appears to be converging towards zero, which 
was first predicted in Ref . Il5l and can be well understood 
using the plasma analogy [iq . 

Even though the relative statistics of two MR quasi- 
holes are trivial, they do pick up a phase due to their 
wrapping around the electrons, analogous to what oc- 
curs in the Laughlin case. Fig. |5] shows that as the size 
of the system increases, the phase accumulated by inter- 
changing two quasiholes (filled symbols) or braiding one 
around the other (open symbols) can be well approxi- 
mated by assuming the wave function rotates by tt for 
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FIG. 2: For X > (x < 0) filled symbols show the phase 
accumulated by interchanging two quasiholes around a circle 
with opening angle ^ centered on the equator (north pole), for 
various A'^ as in Fig. The straight dashed lines in the upper 
half are 0.5(A'' -I- 1/4)(1 — cos^), corresponding to the expec- 
tation of the number of electrons enclosed by the loop. The 
-1-1/4 accounts for the charge pushed out by one of the quasi- 
holes. For X < 0, the dashed lines are — 0.5(A'^ -I- 1/4) cos^, 
i.e., one half of the number of electrons inside minus one half 
the number outside the loop. Open symbols, corresponding 
to a similar calculation with one quasihole moving and the 
other fixed at the center of the circle, almost overlay the filled 
symbols, confirming the trivial relative statistics. 

each enclosed electron (compare to 27r for the Laughlin 
state), when the poles are not separated by the loop (and 
the effect of the pole singularities is analogous to that in 
the Laughlin state). Even for systems consisting of only 4 
electrons, this approximation stays quite good if we cor- 
rect the average electron density for the charge pushed 
out by one localized quasihole (see dashed fines in Fig.[2J). 
This method of correcting the average density also works 
for the Laughlin state on the sphere. 

We now turn to 2n = 4 MR quasiholes, which is the 
simplest case when statistics can be non-Abelian (the 
ground state has degeneracy 2). While the above results 
for 2 quasiholes are anticipated by the plasma analogy 
[lij , one may need deeper CFT jj, 0| arguments in or- 
der to understand the following findings. In the calcula- 
tion, we first fix all quasiholes on the equator and then 
interchange an adjacent pair of them around a circle with 
different opening angles ^ centered on the equator. Pa- 
rameterizing a unitary matrix U by 

f e*''cos/3/2 ie~''/^smf3/2\ 
\ie"'/^sm(3/2 e-"'cos/3/2 )' ^ ' 

we plot in Figs.|21and^the results (in a convenient basis) 
for the transformation Ui corresponding to the braiding 
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FIG. 3: Same as the upper half of Fig. |5| but now with four 
quasiholes present, two of which are fixed on the equator, 
at <j!> = ±37r/4, and two interchanged, with initial and final 
positions ai (j) = ±^ on the equator. The straight dashed lines 
are 0.5(N + 3/4) cos^ - 1/4. Here, +3/4 accounts for the 
average electron-density correction for the charge localized at 
2n — 1 quasiholes. The additional phase offset of —1/4 refiects 
the Abelian part of the braiding statistics, in agreement with 
the predictions of Refs. [Mll3||. 



operation on one of the quasihole pair. Due to the ro- 
tational symmetry around the vertical axis, knowing Ui 
we can deduce other transformations U2-, C^a, and U4 (for 
interchanges of pairs ordered along the equator) by rotat- 
ing and projecting the initial basis and correspondingly 
transforming Ui. It is then easy to show that Ui = U3 
and U2 — U4 due to the form ^ of the wave function. 
Furthermore, we find numerically that U2 ~ F'^UiF, 
where F — {az — ax)/V^, cr's being the usual Pauli ma- 
trices. This approximation is good within a few percent 
for smaller systems and is even better for larger ones. 

According to Fig.01 we see that apart from the Abelian 
phase x^ Ui can be approximated by Ui ~ diag(l + i, 1 — 
i)/\/2, with the disagreement becoming smaller for larger 
systems. Using F, we can then construct all other ma- 
trices Ui- After performing the above approximations. 




FIG. 4: Parameters -q and /3 defining transformation matrix 
@ for the same operations as x shown in Fig. |3 The dashed 
line shows 1/4, an approximation used for in the text. Sim- 
ilarly /3 can be approximated as zero [so that e in Eq. Q is 
not defined]. These approximations become better with larger 
system size and for intermediate cos ^ when the quasiholes re- 
main further apart. The symbol convention is the same as in 
Fig. El Lines interpolate Monte Carlo results. 



we find that the unitary transformations corresponding 
to the braid operators realize the right-handed spinor 
representation of S0(2n)xU(l) (restricted to 7r/2 rota- 
tions around the axes) as predicted in Ref. [13 using 
CFT. In addition to the usual relations required of a 
representation of the braid group on the plane, on the 
sphere the generators must obey an additional relation. 
For the case of 2n = 4, for example, we expect to have 
U1U2U3U3U2U1 = 1. One can easily show that (for gen- 
eral n) the relevant representation of the braid group pre- 
dicted in Ref. satisfies this additional relationship up 
to an Abelian phase. (The failure of the Abelian phase 
to satisfy this law is related to the gauge singularities, 
and will be discussed elsewhere.) 

In summary, we formulated a numerical method to 
study braiding statistics of FQH excitations and applied 
it to perform the first direct calculation of the non- 
Abelian statistics in the MR state. Our findings confirm 
results previously drawn within the CFT framework. 

We have enjoyed helpful discussions with B. I. 
Halperin, C. Nayak, F. von Oppen, and N. Read. This 
work was supported in part by NSF Grant DMR 99- 
81283. The computations were performed on the O2000 
supercomputer CPU farm at Boston University. 
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